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LINEAR RECURRENCE SEQUENCES WITH INDICES IN 
ARITHMETIC PROGRESSION AND THEIR SUMS 

DANIEL BIRMAJER, JUAN B. GIL, AND MICHAEL D. WEINER 


Abstract. For an arbitrary homogeneous linear recurrence sequence of order d with con¬ 
stant coefficients, we derive recurrence relations for all subsequences with indices in arith¬ 
metic progression. The coefficients of these recurrences are given explicitly in terms of 
partial Bell polynomials that depend on at most d — 1 terms of the generalized Lucas 
sequence associated with the given recurrence. We also provide an elegant formula for the 
partial sums of such sequences and illustrate all of our results with examples of various 
orders, including common generalizations of the Fibonacci numbers. 


1. Introduction 


Let d be a positive integer and let ( a n ) be a sequence satisfying the recurrence relation 

a n = cia n _i H-b c d a n _ d for n > d, c d / 0. (1.1) 

While it is not surprising that any subsequence of the form (a mn + r ) ne N, for fixed m € N 
and r £ No = N U {0}, also satisfies a linear recurrence relation of order d, little is actually 
known about the structure of the coefficients of these recurrences. In this paper, we answer 
this question in full generality and give explicit formulas in terms of partial Bell polynomials 
in the coefficients ci,..., c d of the original recurrence relation. 

To this end, we introduce the associated sequence 

do = d, a n = V'---{ T R nfe (l!ci,2!c 2 ,...,d!c d ,0,...) for n > 1, (1.2) 

■' (to — 1 ! 

k= 1 v ’ 

where B n ^ = B n ^(x i, x 2 ,...) denotes the (n, /c)-th partial Bell polynomial in the variables 
xi,X 2 , ■ ■ ■ ,x n _fc + i. These polynomials, introduced by Bell [I], provide an efficient tool to 
work with linear recurrence sequences and their convolutions. For the definition and basic 
properties, see e.g. 13 Section 3.3]. 

It can be shown that ( a n ) satisfies the same recurrence relation as ( a n ). For the special 
case of the Fibonacci sequence ( F n ), where d = 2 and c\ = c 2 = 1, the associated sequence 
(F n ) is given by 


Fq = 2 , 


n 


fc=l 


(fc-1)! 

(to — 1)! 


£„,*(!!, 2 !, 0 ,... 



This is precisely the Lucas sequence [9] A000032]. Moreover, if (a n ) is the generalized 
Fibonacci sequence of order d (with c\ = c 2 = ■ ■ ■ = c d = 1), then (a n ) is the corresponding 
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generalized Lucas sequence studied in [6], For this reason, we call the sequence defined by 
(11.21) the Lucas transform of (ci,..., c<j). One of the main features of a n is that it can be 
written as 

d 

a n = ^2 for n > 0, (1.3) 

3 = 1 

where the ctj’s are such that (1 — ait) •••(! — a d t) = 1 — cit — • • • — Cdt d . The equivalence 
of (11.21) and (11.31) was observed by the authors in pj. 

The main result of this paper (see Theorem 12.1|) is that for an arbitrary linear recurrence 
sequence with constant coefficients ci,..., c f i, as given in (11.11) . and for any fixed m £ N 
and r € No, the subsequence (a mn ^r) n ^ satisfies the linear recurrence relation 

a mn+r Tl ®m(n—l)+r T 72 ®m(n-2)+r T ' ' ' T 'Yd Q , m(n—d)+r for n — d, 

with 7 fc = Y!j =i B k:j (0\a m , 1 \d 2m ,..., (fc - j) ! «(fe-j+i)m) for fc = 1,..., d, where (a„) 

is the Lucas transform of (ci,..., Cd)- 

In Section [2j we will prove this result and will illustrate our formula with examples of 
recurrences of order 2 and 3. We will also consider convolved Fibonacci sequences whose 
characteristic polynomials have roots of higher multiplicity. For brevity in our exposition, 
the number of examples discussed in this section is rather limited. However, all of the results 
presented in this paper are valid for homogeneous linear recurrence sequences of arbitrary 
order with constant coefficients over any integral domain. 

In Section [3l we turn our attention to the partial sums of a general linear recurrence 
sequence (a n ) with characteristic polynomial q(t) = 1 — c\t — • • • — c d t d , and give an elegant 
formula for J2J=o a j t erms of a n +i, • ■ ■, a n +d, see Theorem 13.61 To this end, we first 
consider the sequence (y n ) with generating function 1 / q(t) and find a formula for its partial 
sums. The sequence (y n ) is the invert transform of (ci,..., Cd), and together with the 
sequences with generating functions t 3 /q{t) for j = 1 ,d — 1, they generate a basis for 
the space of linear recurrence sequences of order d with coefficients ci,..., Cd, cf. [2] or [10] . 
The formula provided in Theorem 13.61 is carried out for several basic examples. 

Because of the explicit nature of our two theorems, they can be easily combined to find 
formulas for sums of the form )T)” =0 a m j +r . This is discussed at the end of Section [3] for 
recurrence sequences of order 2 and 3. For illustration purposes, we finish the paper with a 
few examples concerning the Tribonacci sequence. 

2. Indices in arithmetic progression 

Let (a n ) be a sequence satisfying the recurrence relation dm and let (a n ) be the Lucas 
transform of the coefficients (ci,..., c ( i), as defined in (11.21) . We start this section by showing 
that a n admits the representation (11.31) . Let a \,..., ad be defined by (1 —ait) • • • (1 — a^t) = 
1 — cjt — ■ ■ ■ — Cdt d , and let s n = i a j f° r n — 0- 

In [21 Proposition 7], the authors showed that for n > 1, 

= V(—ir +fe ^^H nife (l! ei , 2!e 2 ,..., d!e d , 0,...), 

ti 
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where ei,..., e<j are the elementary symmetric functions in 07 ,... , oq. Since ej = (—1 y +1 Cj 
for every j = 1,,d, the homogeneity properties of the partial Bell polynomials give 

B n ,k{ l!ei, 2!e2,..., died, 0 ,...) = (—l) n+fc B n fe (l!ci, 2!c2,..., died, 0 ,...), 

which implies 


Sn 


£ y——fc(l!ci, 2!c 2 , • • •, d!c rf , 0 ,...) 
rl (n- !)! 


a n for n > 1 . 


Since so = ao, we conclude that s n = d n for all n, as stated in the introduction. Using the 
representation m, it is clear that (d n ) satisfies the same recurrence relation as (a n ). 


Theorem 2.1. Let (a n ) be a linear recurrence sequence of order d > 1, satisfying the 
relation a n = cia n -i + • • • + Cda n _d f or n- > d, Cd ^ 0. Let ( a n ) be the Lucas transform 
of (ci,... ,Cd)- For any fixed m € N and r € No, the subsequence ( a mn+r ) ne pj satisfies the 
linear recurrence relation 


Q”mn+r 7l ®m(n—l)+r “1“ 72 ®m(n— 2 )+r T ' ' ' T 'fd &m(n—d)+r ™ — d, 

where each 7 *, is given by 

Ik = / y Bk,j(0ld m , ■ ■ ■ ,{k j)ld(k—j+l)m)■ 

3 =1 

Proof. Since the sequences (o m n+r)neN and (a mn )ne n satisfy the same recurrence relation, 
it suffices to consider the latter. Using the representation (11.31) . for m £ N, we get 


a mn = £ a T n = £k- 

j= 1 3= 1 


thus for n > d, (a mn )neN satisfies the recurrence relation 

- _ (m) ~ _ (m)~ . . /■_1 ( m ) ~ 

®mn a m (ri— 1) ®m(n— 2 ) i i ( -LJ & m(n—d)i 

where ..., e£, are the elementary symmetric functions in a™,..., a™. 

For every k = 1,..., d, let 7 *, = (— l) k+1 e^. Once again, by [21 Proposition 7], we have 


O'm.n , — 


k =1 

n 


£ _ ^' ! - g n,fc( l! 7i ) 2! 72i • ■ •, 0,...), 

fc=i 


(n- l)!a mn = £(fc - l)!S nifc (l! 7 i, 2 ! 7 2 ,...,d! 7 d, 0 ,...). 

fc=i 


and therefore 
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Finally, Lagrange inversion gives 

^ (-i)J-t- 1 

Ik = / y F k j(0 !®mi ■ ■ ■ , (k j) !o(fc_ j+l)m) • 


i=i 


This proves the claimed recurrence relation for the sequence ( d mn ) ne ^, and therefore for 
any sequence of the form (amn+rjngfsj. □ 


Remark. Clearly, 71 = d m , and since 7 = (—l) d+i ef l> , we have 7 ^ = (—l)( d+ 1 )( m+ 1 )c™. 
Here is a basic example: 

Example 2.2 (fc-Fibonacci). For /c G N let {F kj7l ) n& ^ be the sequence defined by 

Fk ,0 = 0, F k> 1 = 1, and F fcin+ i = kF k ^ n + F k . n _ x for n > 1. 

In this case, (F k)n ) ne ^ is the fe-Lucas sequence denoted by (Lfc, n )neN in the existing literature 
(see e.g. 0 ). By means of Theorem 12.11 we then get 

Fk,mn+r ^k,m F kj1n (n—l)+r T ( 1) 2)+r ^ — 2. (2-3) 

Moreover, the representation (II.2p gives the identity 


(J - 1)! 
^ (m- 1)! 


m—1 


E 

i=0 


m I m — j 


F k ,m — } ^ . \ - B m ,j (l!^, 2!, 0,. 

• I m — J_) J ^^ jji — j \ j 


k m ~ 2 i 


The recurrence relation (|2.31) coincides with the one given in [5j Lemma 3]. It is easy to 
check that L k ,m — 1 T 1 * 


Remark. An interesting consequence of Theorem 1 2. II is that the structure of the recurrence 
relation satisfied by any arithmetic subsequence of a given linear recurrence sequence with 
constant coefficients, only depends on the order of the given recurrence. For example, for 
any linear recurrence sequence (a n ) of order 2 with coefficients ci,C 2 , we always have 


O’mn+r — ®m(n—l)+r T ( 1) C 2 Om(n— 2 )+r ^ 2, 

and for a linear recurrence of order 3 with coefficients ci, C 2 , C 3 , we get 

®mn+r O m ( n —i)_|_ r T 2 (® 2 ?n ® m ) CL m (n-2)-t-r T C 3 3)+r for 71 A 3, 


where (a n ) is the Lucas transform of the coefficients of (a n ). Thus the key is to understand 
the terms a m , CL 2 mi ■ ■ ■, a*n) for which the representation in terms of partial Bell polynomials 
may be useful. 


In order to illustrate the use of (11.21) , we now consider two examples of linear recurrence 
sequences of order three. They both use the following identity: 


B 


m,j 


(*„ * 2 . *s, o, ■ ■ •) = E t (A) LSE) (»)'(§) 


3j—m—2i 


e=o 


fX 3 _\m- 2 j+e 
\ 3! 1 
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Example 2.4 (Tribonacci, A000073 in |9j). Let (t n ) be defined by 

to = h = 0 , t2 = 1 , 
tn = t n —i + t n —2 + t n —3 for n > 3. 

Theorem IQ gives the recurrence relation 

tmn+r = t m t m ( n - i)_|_ r + ~ 2)+r + ^m(n—3)+r f° r n — 3, 

where (t m ) is the Lucas transform of (1,1,1). By (11.21) . we have 

m , . 1 . m -1 .i 


(2.5) 


tm — 


^ (m — 1 )! 

7 = 1 


m — j 


nn ( m — j 
i 


3-t 


J=1 ' ' 3=0 i=\j/ 2] 

This is sequence [9j A001644] and can also be described by 

to = 3, ti = 1, 4 = 3, and t n = t n -1 + tn- 2 + 4-3 for n > 3. 

The recurrence relation (|2.5I) is consistent with the one obtained in [7] Theorem 1], 
Example 2.6 (Padovan, A000931 in [9]). Consider the sequence defined by 

Po = 1 , Pi = P2 = 0 , 

P n = Pn -2 + Pn -3 for 71 > 3. 

Theorem m gives the recurrence relation 

Pmn+r Pm Pm{n— l)+r P 2 (P^m Pm)Pm(n— 2)+r T P m {n— 3)+r ^ — 3, (2-7) 

where (P n ) is the Perrin sequence [9, A001608]. It satisfies the same recurrence relation as 
(P n ) but with initial values Po = 3, P\ = 0, and P -2 = 2. Moreover, by (11.21) . we have 

m / . . x, m— 1 

U - !)• 


p — 

4 rn. — 


(m — 1 )! 


Pmj(0,2!, 3!, 0,...) = ^ 


m 


m — j 
m — j \2j — m 


(2.9) 


J=1 ' ' j= Tm/21 

Example 2.8 (Narayana’s cows sequence, A000930 in [9]). Let (N n ) be defined by 

Nq = Ni = N 2 = 1 , 

Nn = An -1 + An -3 for n > 3. 

Once again, by Theorem 12.11 we get the recurrence relation 

N mn +r — Nm A m ( n _ 1 ) +r + 2 (A 2 ?n A?n.) A m ( n _ 2 )+r T A m ( n _ 3 )_|_ r for 71 > 3, 
where 

= v -EAbwi. 0,31.0,...) = l( ”y; /2J 

^ m-1! JV ’ ’ ’ ’ ' ^ m — 2j \ j 

3=1 K J 3=0 J V J 

While (A n ) counts the number of compositions of n into parts 1 and 3, it can be shown that 

(A 3 n _i) counts the number of (” 2 X )-color compositions of n. Since A 3 = 4 and Nq = 10, 

this subsequence satisfies the relation 

A 2 = 1 , A 5 = 4, Ns = 13, 

A 3 n +2 = 4A 3 ( n _ 1 ) +2 — 3 A 3 (n _ 2) +2 + A 3(n _ 3)+ 2 for n > 3. 
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We finish this section with a linear recurrence sequence of order 4 whose generating 
function has roots of multiplicity 2. 

Example 2.10 (Convolved Fibonacci, A001629 in [9j) . Let (a n ) be the sequence obtained 
by convolving the Fibonacci sequence with itself. This sequence can be described by 

oo = a i = 0, «2 = 1, a,3 = 2, 

0, n — ^^n—l + Qn—2 2tt n—3 ^n —4 for U ^ 4. 

In this case, the Lucas transform a n of (2,1, —2, —1) satisfies a n = 2 L n , where (L n ) is the 
Lucas sequence [9J A000032]. By Theorem 12.11 for n > 4 we then get 

7l®m(n—1) 4“ l2&m(n— 2) T T3®m(n—3) 4“ 74®m(n—4) 

with 

71 — — 2L m , 74 — 1, 

72 = 2 (^ 2 ™, — h m ) = L / 2 m ~ 2L m = 2(—1) + — L m , 

73 = g(2fl3m 3u m Q.2m. T O m ) = 3 (L 3m — 3L m Z/2m 4~ 2 L m ) = ( 1) 2 L m . 

Here we have used the known identities L 2 m = — 2(—l) m and L 3m = — 3(—l) m L m . 

In conclusion, for n > 4 we have 

O’mn 2 L m CLm^n—l) (2( 1) 4“ 2) 4~ ( 4) 2-L m fl m ( n _ 3 ) ®m(n—4) - (2-11) 

For the special cases m = 2, 3,4, 5, we have L 2 = 3, L 3 = 4, L 4 = 7, L 5 = 11, and so 

«2n = 6 a 2 ( n _i) — Hfl2(n—2) + 6 a 2 ( n _ 3 ) ~ a 2(n-4); 

« 3 n = 8 a 3 ( n _!) — 14a 3 ( n _ 2 ) — 8 a 3 ( ra _ 3 ) — a 3 ( n _ 4 ), 
a 4n = 14a 4 ( n _!) — 51a4(„_ 2 ) 4- 14a 4 ( n _ 3 ) — a 4 ( n _ 4 ), 
a 5 n = 22 a 5 ( n _ 1 ) — 119a 5 ( n _ 2 ) — 22 a 5 ( n _ 3 ) — a 5 ( n _ 4 ). 


3. Sums of linear recurrence sequences 

For fixed ci,..., iq with c d / 0, let 

<?(f) = 1 - cH - c 2 t 2 - c d t d , (3.1) 

and let ( y n ) be the sequence with generating function Y(t) = l/q(t). Denoting Co = —1, we 
then have 

/ d \ / 00 

1 = q(t)Y(t) = ( - Cntn ) ( VntT 

' n =0 ’ ' n=0 

which implies c* Vn—i = 0 for every n > 1. Therefore, 

dsn \ d n d , j 

1 — ^ ^ ^ ^ CiUn—i j — ^ ^ ^ ^ CiVn—i ^ ^ ^ ^ Q j Ud—j 

n =1 ' z=0 ' 72—0 z=0 j=0 ^ z=0 

and so 

d \ d —1 / j \ d —1 / d 1 j 

J^cijyo = 1+5^( J^cijyd-j = i+ 5 Z( Ci )yj+ 1- 

2=0 j =0 ' 2=0 j =0 ' 2=0 




(3.2) 
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This is the base case for the following statement. 

Proposition 3.3. Let (y n ) be the linear recurrence sequence with generating function l/q(t), 
where q(t ) = 1 — c\t — C 2 t 2 — • • • — Cdt d with Cd ^ 0, and let Co = —1. Then for n > 0, 

n d— 1 / d— 1— j \ 

= 1 + J2 ( J2 c i\lJn+j+l■ (3-4) 

j =0 j =0 ' j=0 

Proof. We proceed by induction on n. The base case n = 0 was established in (13.211 . Assume 
that (13.411 holds for n — 1. Then 

n n —1 d —1 , d—l—j 

Vi = 9(!) y i + ^( 1 )^ = 1 + £ £ CiUn+j + q(l)Vn 

j =0 j =0 j =0 ' i=0 

d ~l ,d-l-j s d -2 , d—2—j . 

= 1 T ^ ^ ' Cj j Un+j — C-dVri = 1 T ^ ^ ^ ) Cj j T/riH-jH-1 CdUn 

7=1 A i=0 ' j=0 A i=0 ' 

rf _2 , d—l—j . d- 2 

= f T ^ ] f 'y ] Cj j ^ ' C,7—1— 7 'yn+ 7 +l — c d2/n 

7=0 A i=0 ' 7=0 

d—2 ,d—l—j s d —1 ,d—l—j n 

= 1 + ^ ] f ^ ] c « J 2 /n+ 7 +l — Un+d = 1 + ^ ^ ^ ] c i J Un+j+l- 

7=0 A i =0 ' 7=0 A j=0 

Hence the identity (13.411 holds for all n > 0. □ 

Let q(t) be as in (13.111 . For l € {0,1,..., d — 1} we let (yiP) be the linear recurrence 

sequence with generating function Yg(t) = t l /q(t). Note that (yPp is the sequence (y n ) 
introduced above, and for f>0we have 

y { 0 ] = ■■■ = y{-\ = 0 and Vn' ] = Vn-t for n > i. 

Clearly, the sequences (yPP ), (ylP), ..., (ylf form a basis for the space of all linear 
recurrence sequences of order d with coefficients c\,... ,Cd- 

More precisely, if (a n ) is a linear recurrence sequence satisfying a n = cia n _i + - • • + Cda n -d 
with initial values ao,..., ad- 1 , then 

a« = A 0 yi 0) H-1- A d _iy^ _1) , where 


Ao = ao and \ n = a n — c j a n-j for n = 1,. .., d — 1. 

i= 1 

Theorem 3.6. Let (a n ) be a linear recurrence sequence of order d satisfying 

a n — a.\a n —\ T * * * T Cda n —d for n P d, 
with initial values ao,..., ad- 1 , and let Co = —1. For n > 0, we have 

n d —1 /d—l—j v 

?(f) J2 a i = J2 ( Ci ) { a n+j+i - a P >. 


(3.5) 


7=0 


7=0 v i=0 


DANIEL BIRMAJER, JUAN B. GIL, AND MICHAEL D. WEINER 


where q{ 1) = 1 — c\ — • • • — c^. 

Proof. We start by writing cij = Aq yf' 1 + • • • + Xd-iyj 1 ^ as in (J3.5I) . Thus 


n d— 1 


d-1 


d-1 


n—t 


d(l)X a l = dWXX^X = q{ J^yj-t) = XM Q^^Vj 

j=o j =o i=o e=o A j =£ J f=0 v j—Q 

which by (13.41) becomes 

n d— 1 / d— 1 ✓ d—l—j \ 

^X^' = XM 1 + X ( X Ci)yn+j+l-£ 

j=o e=o A J=0 v i=0 / 

d—1 d—1 d—1 / d—l—j . 

= X Xi + XX ( X Ci ) 2 / "+j+ 1 -^ 

^=o £=o j=o ' i=0 ' 

d— 1 d—1 , d—l—j v d—1 d — 1 d—1 , d—l—j 

= X A < + X( X C 0 X Al2/n+j+l-l = X ^ + X ( X c ij 0 n+j+l- 

f=o j =o A j=o A £=0 £=0 J=0 \ i=0 

Now, by means of (I3.5|) . we have 

d—1 d—1 j d— 1 d—l—j 

X^ = -X(X%k-w = -X( X c *)°t 


e=o 


j =0 i =0 


j=0 i=0 


and therefore, 


n d— 1 .. d— 1 —j , 

9 ( 1 )X“i = X X Ci )( a «+j+i - a j)> 

j=o j=o A i =o A 

as claimed. 

Example 3.7 (d-step Fibonacci). Let d € N with d > 2. Let (fjfi) be defined by 

A’ = ■ ■ ■ = & = 0. A-i = 1. flf’ = fn -1 + ■ " + & for n>d. 
By Theorem 13.61 


□ 


n d— 1 

Ad) _ i 


E ft = ih E^ - 2 - i) 

j =o j =o 


■(d) 

n+j+1 


d—1 


-/f) = T ij(E(‘ i - 2 -j)h« + i 

A i=o 


+ 1 


Example 3.8 (d-step Lucas). Let d € N with d > 2. Let (4f^) be the L-sequence associated 
with (fn^) ■ It satisfies the recurrence relation 

£$ = d, lf ] = 2 j - 1 for j = 1,..., d - 1, 

4 d) =4- 1 + -"+4-d for n > d. 


By Theorem IX 


d-i 


E4 d> = TbE<‘ i - 2 -j> 

3=0 j =0 


? (d) _ Ad) \ 

'n+j+1 L? 1 ’ 
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which can be written as 


E<i E(^ 2 U)Wi 

3=0 V 3=0 

In particular, for d = 2 and d = 3, we get 


d -1 


(d) d(d-3) 


(3.9) 


E 4 2) = Cj 2 - 1 and £ if = i (« 3 - Cii) = * (4+2 + 4 3) ), 
3 =0 3=0 

which are sequences A001610 and A073728 in [9], and for d = 4, 

n 

i + 4?’ + 2 ). 

3=0 


Subsequences with indices in arithmetic progression. As discussed in Theorem l2.ll 

given a linear recurrence sequence (a n ) with constant coefficients, any subsequence of the 
form (a m n+r)neN also satisfies a linear recurrence relation with constant coefficients that 
depend on (a n ), the Lucas transform of the coefficients of (a n ). Consequently, Theorem 13.61 
may be used to derive, in a straightforward manner, formulas for the sums J2j=o a mj+r- 
In order to illustrate the combined use of these theorems, we will discuss some examples 
for linear recurrences of order two and three. The higher the order of (a n ), the more terms 
of the associated sequence (a n ) are required to find the coefficients of the recurrence relation 
satisfied by (a mn + r ) nS N. However, the number of terms needed is one less than the order. 
More precisely, if the order of (a n ) is d, we will only need to compute a m , a 2 m , ■ ■ ■, &(d-i ) m - 

Example 3.10 (Linear recurrences of order 2). Let (a n ) be defined by 

o, n — c\Q n —i T C 2 d n —2 for 71 ^ 2, 


with initial values ao and ai. By Theorem 12.11 we know 

Q"mn+r ®m(n—l)+r T ( 1) C 2 ®r 7 i(n—2)+r for 71 2, 

where d m is given by 


Clm — 




m— 1 


m ( rn — j\ m _ 


3 =0 


m-j \ j 


(rn-l)\ 

Moreover, by Theorem 13.61 

( ® J m(n-\-2)+r ^ra+r) iflm l)(^m(n+l)+r ^r) 


2 V 0 . 


1 2 


^ ^ timj+r 
3=0 


Om + (- l) m+1 ^-l 

( 1) ^2 ^ran+r T (flm l)d r 

a m - (-l) m c^ - 1 * 


hor the special case ui me 

n 


E j-, -^fc,m(n+l)+r ( 1) ^l,mn+r T 1 

J= „ t - m ' +r = - (-1)” -1 ’ 
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and for the fc-Lucas sequence, we have 

l-'k,m(n+l)+r ( 1) ^i,mn+r + {Lk ,ra 


n 

^ ^ Lk^mj+r 
3=0 


Lk,m-(~ l) m -l 
These formulas are consistent with the ones given in m- 
Example 3.11 (Linear recurrences of order 3). Let (a n ) be defined by 

CLn — C-lC^n-l T ^2^71—2 T Qj&n—3 ^r Tl ^ 3, 
with initial values ao, ai, and < 22 . By Theorem 12.11 we have 

®mn+r ^m(n—l)+r T 2 (^ 2 m ^m) ^m(n— 2 )+r “h C 3 U m ( n _ 3 )_|_ r fol n ^ 3, 

where a m = -B mj (l!ci, 2!c 2 ,3!c 3 ,0,...). Theorem [376] then gives 

n 2 / 2 -i s 

ttmj+r — EE Ci j (fl'm(n+j+l)+r ®mj+r) • 

i=0 j=0 ' t=0 ' 

where c 0 = -1, ci = a m , c 2 = \{a 2m - a^), and q( 1 ) = 1 - a m - l(a 2m - a^) - eg 1 . 
For the special case of the Tribonacci sequence (cf. Example 12.411 

to — 1 1 — 0, t 2 — 1, t n — t n — i T t n —2 T t n —3 for n ^ 3, 
the above formula (13.1211 gives 

^m(n+l)+r T (l + 2^t^ m ~ ^m))^mn+r 4" ^m(n—l)+r T ^m,r 


(3.12) 


^ ^ tmj-\-r 
3=0 


tm T 9 (bm ^m) 


where 7 mir = (t m + ^(t 2m — — l)t r + (im — l)i m +r ~ t 2m+r . Here are a few values of the 

sequences ( t n ) and (t n ), taken from [9]: 

(A000073) tn : 0,0,1,1, 2 ,4, 7,13, 24,44,81,149, 274, 504,927,... 

(A001644) tn : 3,1, 3, 7,11, 21,39, 71,131, 241, 443,815,1499, 2757,... 

Tribonacci numbers have been extensively studied, and some special cases of the above 
formula can be found in the literature, see e.g. [8] and [71 Theorem 3]. 

We finish this section with a short list of particular instances of the above sum. 

n 

^ ^ tj — 2 (tn +2 T t n l) j 

3=0 


t2j — I (t2n+l + t 2n ) , y] t 2 j + 1 — | (f 2n + 2 + t 2n+ \ — l) , 


3=0 

n 


3=0 


^ " t-tj — 2 (^3n+2 ton l) > ^ " t^j — 4 (^4n+2 T ^4n l)j 

3=0 j=0 
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where /q = — 1 . I\ 


^ ' t5j+r — 22 (^5ra+2+r “I - 8t5n+l+r “I - ^>tbn+r “I - Ir) ; 

j =0 

= —9, I 2 = 7, I 3 = —3, and I 4 = —5. 


Acknowledgement. The authors would like to thank James Sellers for bringing sequences 
with indices in arithmetic progression to their attention. 
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